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Abstract This note gives a correction to the standard analysis of the 
delay pattern in the radio signals from a pulsar in a binary system; the same 
coordinate frame should be used for the transmission of the signal as for the 
motion of the pulsar in the field of its companion. 



The binary pulsar has been dubbed a " unique gravitational laboratory" [T] ; 
the title is apt, because we are on the threshold of observing the spin-orbit and 
tidal aspects of the system. But a necessary prerequisite is that we have a correct 
description in the point-particle approximation. The variable delay of a radio 
signal from a pulsar in a binary system as it traverses its orbit has been analysed 
by Damour and Taylor [2], and their results continue to be used in observational 
analysis up to the present time. The purpose of this note is to point out that 
their article contains an error, due to different choices in the coordinate frames 
used in the analyses of the Roemer and Shapiro delays respectively. 

The first analysis of the Shapiro delay was by Blandford and Teukolsky[3] 
(BT), but they did not propose its use in observational analysis, because it 
is of the same order of magnitude as the correction, arising from first post 
Newtonian (lpN) modification of the Kepler orbit, to the Roemer delay. We 
begin by considering the case that the pulsar's mass is negligible compared with 
that of its companion, so that the system in question is planetary. A Kepler 
orbit is specified in the observer's coordinate system by a set of five parameters, 
namely the longitude of the line of nodes fi, the inclination i to the plane of the 
sky, the angle from the line of nodes to the perihelion the semimajor axis a, 
and the eccentricity e. The position of the pulsar on its orbit is specified by the 
polar angle <f>, measured from the line of nodes. Based on an isotropic metric 

ds 2 =fl-^\ dt 2 - (l + ^\ [dr 2 + r 2 (dO 2 + sin 2 4>)} , (1) 

where terms of order m 2 have been discarded, the Roemer delay from a point 
on the orbit is 

Ar = r sin i sin (f> , (2) 

and the Shapiro delay, which takes account of the varying refractive index of 
space due to the companion's gravitational field and the consequent bending of 
the signal's trajectory, was obtained by BT as 



A a = -2m 
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It should be noted that BT also calculated the Einstein delay time Ae, which 
is the varying gravitational red shift, and which is intermediate in order of 



1 



magnitude; As/Ar is of order v 2 /c 2 while Ae/Ar is of order v/c. The total 
delay is the sum Ar + Ae + As, and the maximum value of v/c for pulsars 
observed so far is about 2xl0~ 3 . 

The lpN correction to Ar was given by Epstein[2], the quantity r in ^ 
being given, in the planetary limit, as 



= a[l-(e + Ae)cosx + /cos2x] , (4) 



where 



and 
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Ae = (13 - 2e 2 ) (2 - e 2 ) , /= (13 - 2e 2 ) . (6) 
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While the parameter Ae is just a small change in the ellipticity of the orbit, the 
other parameter / gives a distortion of the orbit. The consequent corrections 
to Ar are indeed, as anticipated by BT, comparable in magnitude with As, 
but because the latter contains Ar/o an analysis of the delay pattern from a 
double pulsar has enabled Kramer and WexpQ to measure a and i separately, 
whereas previous analysis based on expressions for Ar and Ae allowed only 
measurement of asini. 

However, if we use the same isotropic metric for the planetary motion as for 
the signal transmission, it is easy to deduce that to order m the only relativistic 
correction is in the constant precession rate 

(7) 

that is there is no distortion of the elliptic orbit. If, for example, we use the 
harmonic coordinates advocated by Fock[6], as a precise version of such a metric, 
that is 

dg2 = r-m dt2 _ r_+jn dr2 _ 2 + M ^ (g) 

r + m r — m 

then the relativistic orbit in the plane 9 = Tt/2 may be obtained from the well 
known equation for the Schwarzschild radial coordinate R, namely[5] 

dR\ 2 E 2 -l< 2m 




]■> R + J2 R R + 2mR > ( 9 ) 

where E and J are the energy and angular momentum integrals, simply by 
putting R = r + m, giving, to lpN order, 
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Then, denning new constants a and e by 



E = 1 
this factorizes as 
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a 2 (l-e 2 ) 
so, in harmonic coordinates, the orbit is 

a(l- 



(r — a + ae) (a + ae — r) , (12) 



1 + cos ip 



ip = X(<j>-vj) , (13) 



which is the equation of an ellipse precessing, without distortion, at the rate 
((?]). Note that <fi, as always, is measured from the line of nodes, while ip is 
measured from an initial periastron; successive periastrons occur at intervals 
of 27r.When we take account of the slightly different expressions for a and e 
in the Schwarzschild description, it gives an orbit similar to the Epstein or- 
bit. But putting R instead of r into the expression for the Roemer delay gives 
an answer which contains the additional term msinisin^. It may be verified 
that the Shapiro delay, calculated in the Schwarzschild instead of the isotropic 
coordinates of BT contains another additional term which exactly cancels this 
one. The expression given by Epstein[3] and by Damour and Taylor: 2 j5 which 
combines a calculation of the Roemer delay in the coordinates of Einstein, In- 
feld and Hoffman with one of the Shapiro delay in the isotropic coordinates, is 
incorrect. 

In order to obtain the correct delay formula we now extend the above analysis 
to the case of two masses of comparable magnitude. For this system Fock[6| 
derives the orbit from a lagrangian (see his eqn. (81.01) simplified by using 
(81.18)) 

1 n to 3w 4 / 3to* \ v 2 (3m + m*) m 2 m* 9 
L =2 V + 7 + -( 1 -— j + 2r -ft*-2^ (W) 

to give the energy and angular momentum integrals 

E, = \.2 " ■ 3m* \ E X (6m -7m*) 
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These combine to give the orbit equation 



dr 
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where 



Now we define 
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which, putting 
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gives the same orbit equation as in the planetary case, but with r replaced by 
R. Note that a and e reduce to their planetary values on putting m* = 0, and 
E = 1 + Ei. Thus the general r is simply 



a(l-e 2 ) I 
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The rate at which the orbit is described is obtained from the angular mo- 
mentum integral (fTl)|) , leading to 
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Then, substituting the orbit equation for r, 
(l-e 2 ) 3/2 



/ a ' 



mj (l + ecos'i/O 
giving the secular equation 
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where 

X -2tan- 1 ( x l^— -tan^ ) , (26) 

and n is the frequency 

m / _l8m + m*' 



a 3 \ 4a 

This equation may be inverted to give ip, and hence <f> and r, as a function of t. 
Note that both the orbit and its rate of description, in the harmonic coordinates, 
are remarkably similar to their Newtonian forms. We may finally express the 
Roemer delay 

Afl = rsin0sini (28) 

as a function of t, and this is the quantity which must be added to the BT values 
of Ae and As to give the total delay. 
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